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A complete high precision set of deuteron analyzing powers for elastic deuteron-proton (dp) scat- 
tering at 250 MeV/nucleon (MeV/N) has been measured. The new data are presented together 
with data from previous measurements at 70, 100, 135 and 200 MeV/N. They are compared with 
the results of three-nucleon (3N) Faddeev calculations based on modern nucleon-nucleon (NN) po- 
tentials alone or combined with two models of three nucleon forces (3NFs): the Tucson-Melbourne 
99 (TM99) and Urbana IX. At 250 MeV/N large discrepancies between pure NN models and data, 
which are not resolved by including SNFs, were found at cm. backward angles of 6c. m. ^ 120° 
for almost all the deuteron analyzing powers. These discrepancies are quite similar to those found 
for the cross section at the same energy. We found small relativistic effects that cannot resolve the 
discrepancies with the data indicating that other, short-ranged 3NFs are required to obtain a proper 
description of the data. 



PACS numbers: 21.30.-x, 21.45.-v, 24.10.-i, 24.70. -|-s 

A hot topic in present day few-nucleon system studies 
is exploring the properties of three-nucleon forces (3NFs) 
that appear when more than two nucleons {A > 3) inter- 
act. The 3NFs arise naturally in the standard meson ex- 
change picture in which the main ingredient is considered 
to be 27r-exchange between three nucleons along with the 
A-isobar excitation, the mechanism initially proposed a 
half century ago by Fujita and Miyazawa Further en- 
hancements have led to the Tucson-Melbourne (TM) 
and Urbana IX 3NF 0] models. New impetus to study 
3NFs has come from chiral effective field theory (xEFT) 
descriptions of nuclear interactions. In that framework 
consistent two-, three-, and many- nucleon forces are de- 
rived on the same footing The first non-zero con- 
tribution to 3NFs appears in xEFT at the next-to-next- 
to-leading order (N^LO) of the chiral expansion. That 
explains why 3NFs are relatively small compared to NN 
forces (2NFs) and why their effects are easily masked. 
Therefore, it is, in general, hard to find evidence for them. 

The first evidence for a 3NF was found in the three- 
nucleon bound states, "^H and •^He 0, The bind- 
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ing energies of these nuclei are not reproduced by exact 
solutions of three-nucleon Faddeev equations emjaloying 
modern NN forces only, i.e. AV18 6], CD Bonn 7], Ni- 
jmegen I, II and 93 (8J. The underbinding of '^H and 
■^He can be explained by adding a 3NF, mostly based on 
27r-exchange, acting between three nucleons '9l4ll|. The 
importance of SNFs has been further supported by the 
binding energies of light mass nuclei. Ab initio micro- 
scopic calculations of light mass nuclei, such as Green's 
Function Monte Carlo [l2| and no-core shell model calcu- 
lations [13], highlight the necessity of including 3NFs to 
explain the binding energies and low-lying levels of these 
nuclei. 

In addition to the study of 3N bound states, three nu- 
cleon scattering is an attractive candidate for a further, 
more detailed investigation of properties of 3NFs, such 
as their momentum, spin, and/or isospin dependence. In 
elastic nucleon-deuteron (Nd) scattering and in deuteron 
breakup reactions one can measure not only differential 
cross sections but also a rich set of polarization observ- 
ables under various kinematic conditions. The impor- 
tance of 3NFs in Nd elastic scattering was shown for the 
first time in [l3|- Clear signals from 3NFs were found 
around the cross section minimum occurring at cm. an- 
gle 9c. m. ~ 120° for incident energies above 60 MeV/N. 
Since then experimental studies of proton-deuteron (pd) 
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and/or neutron-deuteron (nd) elastic scattering at inter- 
mediate energies have been performed extensively and 
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15-21| and 



have provided precise data for cross sections 
spin observables, such as analyzi ng p owers 
|24| . spin correlation coefficients [24| . and polarization 
transfer coefficients [H, HI] ■ Large discrepancies be- 
tween data and rigorous Faddeev calculations with mod- 
ern NN forces alone have been reported, which are partic- 
ularly significant in the angular region of the cross section 
minimum and at incident nucleon energies above about 
60 MeV/N. In the case of the cross sections, a larger part 
of the discrepancies are removed combining the NN forces 
with a 3NF such as the TM99 ^] or Urbana IX. Those 
results can be taken as clear signatures for 3NF efi'ects 
in Nd elastic scattering and the data themselves form a 
solid basis to test new theories. In contrast to the elastic 
scattering cross section the theoretical predictions with 
3NFs included still have difficulties in reproducing data 
for some spin observables. 

In Refs. [3, [2l| , precise data for the cross section and 
nucleon analyzing powers at 250 MeV/N in pd and nd 
elastic scattering were reported. Contrary to the results 
obtained at lower energies, the discrepancies in the cross 
sections were only partially removed by incorporating 
3NFs and large discrepancies of ~ 40% still remained 
at backward angles 9c.m. ^ 120°. 

The results of Refs. [l^ [2l| present a new challenge to 
be solved and show that an energy dependent study of 
a variety of spin observables should be a good source 
to obtain a consistent understanding of the dynamics 
of nuclear forces in elastic Nd scattering. Especially 
the deuteron tensor analyzing powers should be inter- 
esting, because matrix elements of the 27r-exchange 3NF 
consist of momentum-spin coupled operators. Data for 
the cross section and spin observables in elastic pd/nd 
scattering have been accumulated at incident energies 
< 200 MeV/N. However, only a few precise data ex- 
ist around 200-300 MeV/N. The new accelerator facil- 
ity at RIKEN, RI beam factory (RIBF), where polar- 
ized deuteron beams are available up to 440 MeV/N, 
offers opportunities to make these experimental stud- 
ies. As the first experiment with a vector and tensor 
polarized deuteron beam at RIBF, the measurement of 
elastic dp scattering at 250 MeV/N was performed, pro- 
viding a complete set of data for all deuteron analyzing 
powers iTii, T20, T21 and T22, in a wide angular range 
0^^ = 36°-162°. 

At RIBF the vector and tensor polarized deuteron 
beam was accelerated first by the injector cyclotrons AVF 
and RRC up to 90 MeV/N, and then up to 250 MeV/N by 
the new superconducting cyclotron SRC. The measure- 
ment of elastic dp scattering was performed with the new 
beam line polarimeter BigDpol, installed at the extrac- 
tion beam line of the SRC. A polyethylene (CH2) target 
with a thickness of 330 mg/cm^ was used as a hydrogen 
target. In the BigDpol four pairs of plastic scintillators 
coupled with photo-multiplier tubes were mounted in two 
independent planes, 90° apart in azimuthal angle and 



operated in kinematic coincidence for elastic dp scatter- 
ing. The opening angle of the BigDpol was 8°~70°. The 
solid angles were determined by the proton detectors and 
the angular spread was A6'iab. ±1°. For subtraction of 
events from the proton knock-out reaction on carbon nu- 
clei ^^C{d,dp)^^B a measurement with a graphite target 
was also performed. 

The data were taken with polarized and unpolarized 
beams for different combinations of the incoming polar- 
ization given in terms of the theoretical maximum polar- 
ization values as {pz,Pzz) = (0,0), (1/3,-1), (—2/3,0) 
and (1/3,1). Those polarization modes were changed 
cyclically at intervals of 5 seconds by switching the RF 
transition units of the polarized ion source. The beam 
polarizations were monitored continuously with a beam 
line polarimeter prior to acceleration by the SRC using 
elastic dp scattering at 90 MeV/N. The analyzing pow- 
ers for this reaction were calibrated in the previous mea- 
surement by using the ^'^C{d, a)^^B* [2^] reaction, the 
(0° ) of which is exactly —1 because of parity conser- 
vation |29j . In the present measurement typical values 
of the beam polarizations were 80% of the theoretical 
maximum values. 

The deuteron analyzing powers for elastic dp scattering 
are expressed through the unpolarized (ctq) and polarized 
(cr) cross sections together with the vector and tensor 
polarizations of the incoming deuteron (pz and pzz) as 

a = (Tq<. I + ViiTii (0) pz sin l3 cos (f) 
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T20 {0)pzz (3cos2/3-l) 



T21 (9) Pzz cos /3 sin (3 sin 
^T22 {9) pz z sin^ PC0S2A , 



(1) 



where and <j) are the polar and azimuthal scattering an- 
gles, respectively f3^. The (3 is defined as the angle be- 
tween the spin direction and the beam direction. In this 
experiment the polarization axis of the deuteron beam 
was rotated with a Wien filter system to the direction 
required for the measurement of a particular analyzing 
power (3l| . For the measurement of the analyzing pow- 
ers «Tii, T20, and T22 the polarization axis was normal 
to the horizontal plane. For the T21 measurement the 
spin symmetry axis was rotated in the reaction plane 
and aligned at an angle j3 — 38.0° ± 0.51° to the beam 
direction. 

Figures [Ija) and (b) show spectra of scintillator light 
outputs of one pair of the deuteron and proton detec- 
tors placed at 19.5° and 55°, respectively. The spec- 
tra obtained with the CII2 target together with the 
graphite target are shown. Impinging deuterons and pro- 
tons punched through the detectors. Identification of the 
scattered deuterons and recoil protons was performed by 
cuts for the spectrum with the CII2 target as shown in 
Figs.[lja) and (b). After the selection of events a spec- 
trum of the time-of-flight difference between the deuteron 
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and proton detectors was obtained (see Fig.[T]Jc)). As for 
asymmetry determination of elastic dp scattering the se- 
lected events within the gate shown in Fig. (Tfc) were 
used. Since the contribution from carbon nuclei was 
small in Fig. [TJc) the background subtraction was not 
performed in the analysis. 

We present these new data and compare them with 
theoretical predictions based on different dynamical in- 
put in Fig.[2j The statistical uncertainties, which are less 
than 0.01 for all the deuteron analyzing powers, are also 
shown. In order to clarify the energy dependence, the 
previously measured deuteron analyzing powers at 70, 
100 and 135 MeV/N [IMS, HI and at 200 MeV/N [H 
are also included together with theoretical predictions. 
The deuteron analyzing powers at 70, 100, 135 MeV/N 
have been reported in Refs. [isl - flTl . [25| in the Cartesian 
notation. Here we show them in the spherical descrip- 
tion. The theoretical predictions shown are based on 
modern NN forces alone and on their combinations with 
3NF models [H, [s^l- We used high precision NN po- 
tentials AV18, CD Bonn, and Nijmegen I and II alone 
[light shaded (blue) bands in Fig. [5] or combined them 
with the TM99 3NF [13 with the cutoff A values that 
yield, for a particular NN force and TM99 combination, 
a reproduction of the '^H binding energy [dark shaded 
(red) bands in Fig. [5]. In case of the AV18 NN force, we 
also combined it with the Urbana IX 3NF j3j] (solid lines 
in Fig. In our calculations we neglected the Coulomb 
force acting between two protons. At energies considered 
here its effect on polarization observables is small [28] . 

For the vector analyzing power iTn the predicted 3NF 
effects at 70 MeV/N are small and the data are satis- 
factorily described, independently of whether 3NFs are 
included. At higher energies, however, a discrepancy be- 
tween data and theory based on NN forces only is clearly 
seen for cm. angles 6'c.m. > 80°. At 100 and 135 MeV/N 
adding TM99 or Urbana IX 3NF provides quite a good 
description of the data. That is not the case at 200 and 
250 McV/N, where only a part of the discrepancy is re- 
moved and the full theory is quite far from data in some 
angular regions. For 9c.m. ^ 60° at all energies 3NF ef- 
fects are negligible and data are reproduced by pure NN 
theory. 

For the tensor analyzing power T20, the predicted 3NF 
effects are smaller than for iTu. Again, at 70 MeV/N the 
theory reasonably well reproduces the data and at 100 
and 135 MeV/N adding 3NFs brings the theory closer 
to the data. At 200 and 250 MeV/N, where the pre- 
dicted 3NF effects for T20 are rather moderate, drastic 
discrepancies between theory and data in the angular re- 
gion 80° < 0c.ni. ^ 140° are clearly seen, which are not 
explained by adding standard 3NFs. 

The data for the tensor analyzing power T21 are repro- 
duced quite well by calculations with the 2NFs only for 
all incident energies. The large effects of the TM99 3NF 
predicted for energies > 100 MeV are not supported by 
the T21 data. They differ clearly from the smaller effects 
of the Urbana IX 3NF, which generally leads to a better 
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FIG. 1: Spectra of light outputs of of the scintillation counters 
for the scattered deuterons (a) and the recoil protons (b) . The 
deuteron and proton detectors were placed at 19.5° and 55°, 
respectively. After selecting events in (a) and (b) the spectra 
of time of flight difference between the deuteron and proton 
detectors was obtained (c). 
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FIG. 2: (Color online) The deuteron analyzing powers iTn, T20, T21, and T22 for dp elastic scattering at 70, 100, 135, 200, and 
250 MeV/N. The light shaded (blue) bands contain predictions of modern NN potentials: AV18, CD Bonn, Nijmegen I and 
II. The dark shaded (red) bands result when those potentials are combined with TM99 3NF, properly adjusted to reproduce 
the '^H binding energy. The solid line is the result obtained with the combination AV18+Urbana IX. The pd data are: at 70 
MeV/N (open circles) from Ref ^2^], at 100 MeV/N (open circles) from Ref. at 135 MeV/N (open circles) from Ref. [13] 
and (solid circles) from Ref. [U], at 200 MeV/N (solid circles) from Ref. [Ij, and at 250 MeV/N (open circles) from the present 
study. 



description for that observable. 

For the tensor analyzing power T22, the discrepancies 
between the data and the predictions based on 2NFs only 
become larger in magnitude and expand to the backward 
angles with increasing incident energy. For that observ- 
able the predicted 3NF effects are especially large and 
similar in magnitude for the 2NFs plus TM99 and Ur- 
bana IX models. At 200 MeV/N and less the predictions 
taking into account 3NFs have good agreement to the 
data at the backward angles, however the data in the 
angular region 40° < ^c.m. S 120° are not described by 
any theoretical predictions. At 250 MeV/N the overall 
agreement is improved by taking into account these 3NFs 
except for the very backward angles. 

All the deuteron analyzing powers, with exception of 
T21, reveal at the highest energy 250 MeV/N and around 
cm. angles ^c.m. ^ 120° large discrepancies to theory 
based on NN forces alone, which are not resolved com- 
pletely by the inclusion of the SNFs. Such behavior of 
the deuteron analyzing powers is quite similar to that of 
the cross section and proton/neutron analyzing powers 
at 250 MeV/N found in Refs. [ll,[2l|. 

The energy dependence of the predicted 3NF effects 
and the difference between the theory and the data for 
the deuteron analyzing powers is not always similar to 
that of the cross section and nucleon analyzing power. 



The vector analyzing power iTn and the tensor analyz- 
ing power T20 have features similar to those of the cross 
section and the proton analyzing power A^. However the 
tensor analyzing power and r22 reveal different en- 
ergy dependence from that of other observables. Starting 
from ~ 100 MeV/N large 3NF effects are predicted. For 
T21 they are of different magnitude for TM99 and Urbana 
IX and the T21 data seem to prefer the smaller effects of 
Urbana IX. For r22 the large effects of TM99 and Urbana 
IX are practically the same. At 200 MeV/N and below 
adding 3NFs worsens the description of data in a large 
angular region. It is contrary to what happens at the 
highest energy 250 MeV/N, where large 3NF effects are 
supported by the T22 data in a large angular range. 

The results obtained for the highest energy of 250 
MeV/N indicate that some significant components are 
missing in the calculations, especially in the regions of 
higher momentum transfer. One possible candidate is 
relativistic effects. We estimated their magnitude for the 
deuteron tensor analyzing powers by comparing nonrela- 
tivistic and relativistic predictions based on the CD Bonn 
potential [s^, [35| . They turned out to be small and only 
slightly alter the deuteron analyzing powers. 

Due to the smallness of the considered relativistic ef- 
fects, it appears that important parts of the SNFs are 
missing. In the meson exchange picture used here, con- 
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tributions of heavy mesons, e.g. ir-p and p-p exchanges 
are omitted. The importance of such contributions is 
also seen in ^EFT, where the leading non- vanishing 3NF 
consists of a 27r-exchange, a Itt exchange-contact and a 
pure contact interaction topology. With increasing en- 
ergy orders of the chiral expansion higher than N^LO 
become important, which introduce a multitude of ad- 
ditional short-range 3NF contributions, coming with dif- 
ferent momentum-spin dependence. The complicated en- 
ergy dependence of the deuteron analyzing powers and 
their strong dependence on the 3NF model used indi- 
cates the importance of such short-range components in 
their full description. The results of Faddeev calcula- 
tions based on a force model with explicit A degrees of 
freedom show large changes in the predicted analyzing 
powers (36j and also indicate the importance of short- 
range 3NF components in the description of these spin 
observables. 

In summary we have reported the first complete set 
of high precision data for the deuteron analyzing powers 
iTii, T20, T21, and T22, in elastic dp scattering at 250 
MeV/N, taken in a wide angular range Oc.m. = 36°~162°. 
These data constitute a solid basis to guide theoretical 
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